Energy transport in strongly disordered superconductors and magnets 
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We develop an analytical theory for quantum phase transitions driven by disorder in magnets and 
superconductors. We study these transitions with a cavity approximation which becomes exact on a 
Bethe lattice with large branching number. We find two different disordered phases, characterized by 
very different relaxation rates, which both exhibit strong inhomogeneities typical of glassy physics. 
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The zero temperature quantum phase transi- 
tions and their quantum critical points have been 
rather well understood in translationally invariant 
systems p. Much less is known about disordered 
systems where the transition is driven by the com- 
petition between strong disorder and interactions. 

Motivated by experiments on disordered super- 
conductors 0, we formulate a theoretical model of 
disorder driven transitions and solve it in the sim- 
plest controlled approximation. Our main results 
can also be relevant for other disordered quantum 
problems, especially for strongly disordered mag- 
nets. The new physics introduced by the strong dis- 
order is the appearance of new phases [S] in which all 
or some excitations are localized in space and have 
infinite lifetime and thus cannot contribute to any 
transport. The quantum critical point at which the 
long range order appears has many features that dis- 
tinguish it from conventional quantum critical points 
in translationaly invariant systems. Most notably, it 
is characterized by a wide distribution of the order 
parameter and the appearance of a new intermedi- 
ate phase in which only low energy local excitations 
have infinitely long lifetime while high energy ones 
can decay. 

Strongly disordered films of InO, TiN or Be dis- 
play a transition from the superconductor to in- 
sulator when their resistivity in the normal state 
exceeds a value of the order of resistance quan- 
tum Rq = 6.5fcfi0. Close to the transition the 
superconductor-insulator phase can also be induced 
by magnetic field; this transition displays a quan- 
tum critical point behavior Q . In the vicinity of the 
quantum critical point, the tunneling spectroscopy 
shows a well defined gap at all points. However, 
the coherence peaks expected for a BCS supercon- 
ductor appear at some locations and disappear at 
others, a phenomenon which is similar to some ex- 
periments in high Tq oxides. The absence of a coher- 
ence peak combined with the detection of a super- 
conducting gap in a single electron tunneling exper- 



iment implies that the disorder does not affect the 
local Cooper pairing of electrons but prevents the 
formation of a coherent state of these pairs. This 
allows to exclude 'fermionic' mechanisms of the su- 
perconductivity suppression through a reduction of 
phonon attraction by Coulomb interactions. Be- 
cause Coulomb interaction is strongly suppressed in 
the insulating phase Q, the most plausible mecha- 
nism for the superconductor-insulator transition in 
homogeneous disordered films is a competition be- 
tween pair hopping and random pair energies on dif- 
ferent sites, as suggested in a seminal paper of Ma 
and LeeQ. 

As we show in this paper the solution of this 
model, which requires going beyond the simple 
mean-field anlaysis of the earlier works, reproduces 
correctly the most important features of the data 
on disordered films: direct superconductor-insulator 
transition, activated behavior close to the quantum 
critical point in the insulating phase, strong depen- 
dence of the activation energy near the quantum 
critical point and huge order parameter variations 
from site to site in the superconducting phase. 

In the absence of Coulomb repulsion the electrons 
are paired even on localized single electron states, 
and pairs can hop from one site to another. This 
physics is described by a Hamiltonian of disordered 
bosons with strong on-site repulsion 0, H, 0] : 

^ = - ( E ^^''^ + E ^^^^ (''t'^j + '^rO I (1) 

Here the state with erf = ±1 corresponds to a lo- 
cal level occupied or unoccupied by a Cooper pair; 
^j's are occupation energies for each site, which are 
quenched random variables drawn from a probability 
P{£,)- Mij describe the pair hopping amplitude be- 
tween sites i and j. These hopping amplitudes cou- 
ple a typical local level to a large number of neigh- 
bors, Z ^ 1. We shall assume that each site is 
coupled to Z neighbours with Mij =^ g/{Z — 1), and 



for technical simplicity we shall study the slightly 
simplified Hamiltonian 



H 




(2) 



but all our conclusions also hold for the case ([T|). 
With a redefinition of the meaning of occupied and 
empty, one can take the site energies to be all 
positive. We shall assume that the ^ are uniformly 
distributed on the interval [0,1/;^]. The important 
feature of this distribution is that it is constant near 
to ^ = 0; the value of v just sets the scale of energies, 
and we can choose v = 1. In this language of Hamil- 
tonians (|ll2p the superconducting phase is mapped 
onto the phase with spontaneous magnetization in 
the X direction, in the insulating phase spins point 
parallel to z-axis. 

The most obvious approach to study this Hamilto- 
nian is through a simple mean-field (SMF) approach, 
where H is replaced by Hmf = ~ Bcrf) 

and B is determined self-consistently. At tempera- 
ture T — 1/(3, this predicts a phase transition from 
insulator to superconductor at the critical value of 
the hopping 
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As P(0) > 0, ^ at low temperatures. 

While these SMF predictions are correct at Z — 
oo, they are qualitatively wrong at low temperature 
in finite connectivity systems. We now turn to a 
more refined mean field discussion, valid for finite 
Z ^ 1, which is the basis for our main results. We 
use a quantum version of the cavity method [lo| 
which would become exact if the spins were on a 
Bethe lattice of connectivity Z. In this method, one 
studies the properties of a spin j in the cavity graph 
where one of its neighbours has been deleted, as- 
suming that the K = Z — 1 remaining neihbours are 
uncorrelated. The system of spin j and its K neigh- 
bours is thus described by the local Hamiltonian 



j^cav 
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Bk<jl + j^aXk] (4) 



where Bk is the local "cavity" field on spin k due to 
the rest of the spins (in absence of j). By solving the 
problem of Z Ising spins in (jH) , one can compute the 
induced magnetization of j, (cj), which is by defi- 
nition equal to Bj/ + ^1 ■ We thus get a map- 
ping allowing to compute the new cavity field Bj in 
terms of the K fields Bk on the neighbouring spins. 
This mapping induces a self-consistent equation for 



the distribution of the B fields [10| . We have made 
one more approximation which is to study the cavity 
Hamiltonian ^ with a mean field method 11]. This 
gives the explicit mapping: 



B, = 



-T 
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In order to understand this mapping, let us imag- 
ine that we iterate it R times on a Bethe lattice. 
For R finite, when the number of spins is large, the 
corresponding graph is a rooted tree with branching 
factor Z — 1 at each node and depth R. The field 
Bq at the root is a function of the fields on the 
boundary. In order to see whether there is sponta- 
neous ordering, we study the value of Bq in linear 
response to infinitesimal fields Bi — B ^ 1 on the 
boundary spins. This is given by 



g tanh(/3^,^ 



K 



(6) 



where the sum is over all paths going from the root to 
the boundary, and the product HnGP '^'^'^^ sites 
along the path P. The response S is nothing but 
the partition function for a directed polymer (DP) 
on a tree, where the energy of each site is e~^" = 
(g/if )(tanh(/3^^)/^„) and the temperature has been 
set equal to one. The solution of this problem, found 
in 'l^ , can be expressed in terms of the convex 

function f{x) = (l/a;)log \k d( {taiih{(3S,) / 

which is minimal at a value x — Xc- In the large R 
limit, there exist two phases for the DP problem: 

• "Self-averaging" (SA) phase: li Xc > 1, then 
(l/i?)logS = /(I) + \og{g/K). The ordered 
phase appears at gc — Ke~^''^^ — g^^'^^ ■ 

• "Glassy phase" (G) phase: If a;c < 1, then 
(l/i?)logS = f{xc)+ log{g/K). The ordered 
phase appears at g^ — Ke"^^^"^^ > 

These two regimes of the DP problem are quali- 
tatively very different. The "SA" regime is the high 
"temperature" phase of the polymer, where the mea- 
sure on paths defined in ([6|) is more or less evenly 
distributed among all paths. The low temperature 
"G" regime is a glass phase where the measure con- 
densates onto a small number of paths. An order 
parameter which distinguishes between these phases 
is the participation ratio Y = X]p'"^pi where is 
the relative weight of path P in the measure ©• In 
the replica formalism the SA phase is replica sym- 
metric, S is self-averaging, and Y = 0; the G phase 
is a one-step replica-symmetry-breaking (RSB) glass 
phase, the value of Y is finite and non self-averaging 



(it depends on the explicit realization of the ^'s even 
in the thermodynamic limit), and its average is given 
by 1 — a^clSl- This glass transition, and the nature 
of the G phase, are identical to the ones found in the 
random energy model 15. lij. 



co*(g)' 




<av>?=0 



g* " g 



FIG. 1: Phase diagram of the spin system for K = 2. 
The right pane shows the critical line Tc{g) while the left 
one shows the critical energy that separates the states 
with zero width from those with a finite width. The 
inset shows the low temperature/energy region. 

Using these DP results one gets the phase dia- 
gram of the spin systems shown in the right pane 
of Fig[T] At any temperature, there is a non-zero 
critical value of the coupling, gc{T), separating an 
"ordered" , superconducting phase with spontaneous 
X magnetization at g > gc(T) from a normal "disor- 
dered" , insulating phase with zero magnetization at 
g < gc{T). Within each phase, there are two regimes 
of temperature, "Self-averaging" and "Glassy" . As 
is clear from our susceptibility analysis, the glass 
transition of the DP affects the propagation of a 
static perturbation in the spin system. In the disor- 
dered SA phase, the total effect of the perturbation 
decreases when R increases, and propagates evenly: 
the average value of the susceptibility coincides with 
its typical value. In the disordered G phase, the 
total perturbation also decays, but it condenses on 
a finite number of paths. Consequently, the sus- 
ceptibility is non-self averaging, similarly to what is 
found in one dimension [l3|. Rare paths are impor- 
tant in the whole "Griffiths" phase where the sus- 
ceptibility distribution has a power law tail. When 
gSMF <_ g <^ the typical susceptibility is finite 
but the average susceptibility diverges. In the or- 
dered phase, the perturbation propagates to infin- 
ity, again with very different patterns in the "SA" 
and "G" phases. The SMF gets the correct result 
of the SA regimes, but completely misses the low- 
temperature physics of condensed correlation paths. 

The RSB transition also strongly affects the scal- 
ing of the field in the ordered phase, for g = gc{l+i)- 
One can study the distribution of fields P{B) in- 



duced by the mapping ([5]). An expansion of it's 
Laplace transform shows that, in the G phase, P{B) 
decays at large B like P{B) ~ C/Bi+^<=. This dis- 
tribution has a diverging mean, dominated by rare 
fiuctuations. A careful analysis of the self-consistent 
equation for P{B) then indicates that the typical 
scale of the field behaves as 



(7) 



In the disordered phase the average value of the 
transverse field is zero, but its quantum fluctuations 
become important. Their main physical effect is the 
broadening of the local levels that, in the absence of 
g, correspond to erf = ±1. At T = the level broad- 
ening means that local excitations of energy w ~ 2^^ 
decay. If we neglect the phonons, energy conserva- 
tion requires that some group of spins with the same 
energy be flipped. In finite systems this is generically 
impossible. To study whether energy can be trans- 
ported, and ergodicity can be restored, in infinite 
systems, we adopt an approach similar to the one 
developed above. Namely, we consider a Bethe lat- 
tice that is very weakly coupled to the environment 
at its boundary and study the effective level width 
of a 'root' spin at a distance R from the boundary in 
the R oo limit. Thus we add to the Hamiltonian 
^ the boundary term Hgnv = — <7jXj{t) where 
Xj (t) are dynamical fields generated by the environ- 
ment, characterized by a response function G{uj). In 
the leading order in g/K the relaxation rate of the 
root spin follows from the Fermi golden rule: 



FoM^ImGH^ n 
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This perturbation-theory-based equation is only 
valid when all fractions inside the product remain 
small 18[ , and the relaxation rate of each spin is very 
small. Thus it is self-consistent if Fq ^ when 
i? — > oo. The typical value of Fo(a;) is controlled by 
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it decreases away from the boundary if fri^^) + 
2\n{g/K) < 0. /r can be computed again by an- 
alyzing a DP problem, which turns out to be al- 
ways in its G phase. At = 0, one finds that 
/r(0)-|-2 ln{g/K) < in the whole insulating regime, 
and the equality is reached at the critical point 
g — gc- One also finds that the region of small — 
gives negligible contribution to /r(a;), which allows 
one to work with the unregularized expression ([8]). 
At non-zero frequencies fri^^) decreases, it is mini- 
mal at ti> = 1/2 (which corresponds to the center of 



the band in our notations). At g < = ii:e/r(i/2)/2 
the relaxation rate is zero for all states, this is the 
superinsulator regime of Q. In the intermediate 
regime 5c (0) > g > g* the states in the middle of the 
band have finite width, they are separated from the 
zero- width states by a critical energy uj*{g) similar 
to the mobility edge of the non-interacting problem. 

In order to understand the low temperature prop- 
erties of the system, it is important to know the de- 
pendence of r(a;) at w > uj*{g). In this regime ^ 
must be modified. Using the mapping to fermions as 
in [isl , one finds the iterative equation for the width 
of the levels on the Bethe lattice: 

r.M^(2,/xf^ J^^j (9) 

This equation is similar to the equation ^ for 
the fields in the ordered phase and can be analyzed 
with the same method, giving the fast level-width 
dependence slightly above co*{g)- 

Ttypiuj) ~ r*e-"°(^')/("-"'(f» . (10) 

This behaviour has important consequencies for 
the low temperature properties of the relaxation, as 
we now discuss. A non-zero but low temperature 
affects the relaxation rate in several ways. First, it 
changes the occupation numbers of the excited and 
ground states, this affects the perturbative equation 
dS]) and thus shifts the position of the ui*{g) line, this 
effect is however small at T <C 1. More importantly, 
a non-zero temperature creates some mobile excita- 
tions with frequencies above u!*{g). These excita- 
tions provide a mechanism for a small level broad- 
ening of the very low frequency levels: They see 
a mobile excitation with energy E with an Arrhe- 
nius rate, giving a width exp{~LJo / {E — uj*) — E /T) . 
The dominant contribution comes from energies E = 
+ V^oT, and results in the temperature de- 
pendence T{g) ~ exp(-2A/cjo/r - uj*{g, T)/T) that 
shows a crossover between a square root and acti- 
vated or even faster behavior as one goes away from 
the critical point. 

In conclusion, we have outlined a solution of the 
strongly disordered spin model on the Bethe lat- 
tice which shows a series of two zero-temperature 
transitions between a phase with no relaxation, a 
phase with a slow relaxation and an ordered phase. 
It also shows that the low temperature phases are 
always very strongly non-uniform: both the order 
parameter formation and the spin relaxation are 
controlled by rare interaction paths containing a 
very small number of spins. When applied to the 
superconductor-insulator transition our results im- 
ply the existence of both weak and strong insula- 
tors. At the critical point the relaxation rate varies 



as exp(l/\/r) but crosses over to activated at lowe4 
g and low T, in the strong insulator the relaxation is 
completely suppressed. Of course, some of the phys- 
ical effects neglected in our model would lead to a 
very slow relaxation even in the strong insulator. 
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